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We present a general method for constructing acoustic analogs of the black hole solutions of two- 
dimensional (2D) dilaton gravity. Because by dimensional reduction every spherically symmetric, 
four-dimensional (4D) black hole admits a 2D description, the method can be also used to construct 
analogue models of 4D black holes. We also show that after fixing the gauge degrees of freedom 
the 2D gravitational dynamics is equivalent to an one-dimensional fluid dynamics. This enables us 
to find a natural definition of mass M, temperature T and entropy S of the acoustic black hole. 
In particular the first principle of thermodynamics dM — TdS becomes a consequence of the fluid 
dynamics equations. We also discuss the general solutions of the fluid dynamics and two particular 
cases, the 2D Anti-de sitter black hole and the 4D Schwarzschild black hole. 
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> ■ I. INTRODUCTION 
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' In recent years it has become increasingly clear that condensed matter systems can be used to mimic various 
kinematical aspects of general relativity 0,12, IS El HI- Analogue models of gravity have been used not only to mimic 
black holes and event horizons but also to describe other spacetime structures, cosmological solutions and field theory 
in curved spacetimes [(J 0, U El El 0, El Condensed matter analogs of gravity could be also used to 

shed light on old and new puzzling features of the gravitational interaction such as the information loss problem in 
the Hawking evaporation of a black hole or the holographic principle. Moreover, analogue models of gravity that 
use condensed matter systems have also the nice feature of being, at least in principle, experimentally testable in 
laboratory. In the near future this could open the way to a "black hole phenomenology" based not on astrophysical 
but condensed matter experiments. 

One strong limitation of this approach is that until now analogue models of black holes can be used to mimic 
the kinematical but not the dynamical aspects of gravitational systems. For instance, one can use acoustic analogs 
to study the kinematics of event horizons, or even the Hawking radiation but not the formation of horizons due to 
\ some distribution of matter or the back-reaction of the geometry on the Hawking radiation (some progress along this 
direction has been achieved in Ref. pdLIT^|) . A further unpleasant consequence of this limitation is the impossibility 
to reproduce black hole thermodynamics using analogue models of black holes. It is well known that the laws of 
bi). black hole thermodynamics are a consequence of the Einstein equations. Although one can define the temperature 
associated with an acoustic horizon 3] , it is extremely problematic to find the right definition of mass and entropy of 
an acoustic black hole. 

At first sight to mimic the dynamical aspects of general relativity using a condensed matter system seems an 
impossible task. The description of the gravitational system is characterized by a huge redundancy of gauge degrees 
of freedom and by a separation between matter sources and gravitational field (respectively described by the stress- 
energy tensor and the curvature tensors in the Einstein equations). Conversely, focusing on acoustic black holes, the 
fluid dynamics is characterized by few physical degrees of freedom (pressure, velocity, density) and does not seem to 
allow for a source-field description . 

Nonetheless, there is a situation in which one could hope to merge gravitational and fluid dynamics. A spherically 
symmetric, black hole with gauge degrees of freedom completely fixed is characterized, owing to Birkhoff's theorem, 
by a handful of physical observables (the charges associated with global symmetries ). These observables could be 
put in correspondence with the fluid parameters. A spherically symmetric black hole can be described by an effective 
two-dimensional(2D) dilaton gravity model |15| obtained from four-dimensional(4D) Einstein gravity by retaining 
only the radial modes. Thus, the most promising framework for trying to merge gravitational and fluid dynamics is 
that of 2D gravity. 

In this paper we will show that this is possible. We propose a general method for constructing acoustic analogs of 
the black hole solutions of two-dimensional (2D) dilaton gravity (Section II). Further, we show that after fixing the 
gauge degrees of freedom the 2D gravitational dynamics is equivalent to a constrained one-dimensional fluid dynamics 
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(Section III). This enables us to find a natural definition of mass M, temperature T and entropy S of the acoustic black 
hole. In particular the first principle of thermodynamics dM — TdS becomes a consequence of the fluid dynamics 
equations (Sections IV and V). The constraint for the fluid dynamics can be enforced in two independent ways. If the 
external parameters for the fluid (profile of the flux tube or potential for the external forces) are given the constraint 
simple determines the barotropic equation of state for the fluid. Alternatively if one chooses a given Equation of 
state the constraint determines the form of the external parameters for the fluid (Section VI). We also discuss two 
particular cases, the 2D Anti-de sitter black hole and the 4D Schwarzschild black hole (Section VII). 



V 



II. 2D ACOUSTIC METRIC AND 2D BLACK HOLES 

Let us consider the generic 2D dilaton gravity model (For a review see Ref. characterized by a dilaton potential 

A = \ f d 2 x^{cj ) R + \ 2 V (<f>)), (1) 



where is a scalar (the dilaton) and A is parameter with dimensions of a length -1 (throughout this paper we will 
use natural units c = fi = ks = 1)- Two-dimensional dilaton gravity models have been used in various situations as 
effective description of 4D gravity . In particular, every spherically symmetric 4D solution can be described in terms 
of a 2D dilaton gravity model. The model JTJ admits black hole solutions, which in the Schwarzschild gauge take the 
form [l6|. 

d,s 2 = - (j(Xr) - dr 2 + (j(Xr) - ™\ 1 dr \ ( f ) = Xr, (2) 

where M is the black hole mass and J = f Vd<f>. The black hole horizon is located at r = r^, with J(Ar^) = 2M/X. 

We want to find acoustic analogs of the 2D black hole solutions |J2J| . Obviously, the most natural candidate is a fluid, 
whose motion is essentially one-dimensional. Let us therefore consider the steady, locally irrotational flow of a 3D 
fluid which is barotropic and inviscid. If the transverse velocities (in the z and y directions) are small with respect to 
the velocity along the x axis, we can set them to zero. The propagation of acoustic disturbances x on the background 
fields is determined, at the linearized level, by the wave equation V 2 % = (1/ ' ^p-g^d^p-gg^ d u x) = 0, where 
is the so-called acoustic metric In the case under consideration, the acoustic metric describes a one-dimensional 
slab geometry where the velocity is always along the x direction and the velocity profile depends only on x, 

ds 2 = ^-[- (c 2 - ug) dt 2 - 2v dxdt + dx 2 ] , (3) 



where pq(x),vq(x) are, respectively, the background density and velocity of the fluid and c{x) — y/ (dP/dpo) is the 
local speed of sound (P is the pressure) . 

The two metrics (|5J|, © can be transformed one into the other using a Painleve Gullstrand-like coordinate trans- 
formation 0, 0] and identifying the function J and the black hole mass M in terms of the physical parameters of 
the fluid. It is convenient to work with a dimensionless fluid density. This can be achieved by multiplying the metric 
@ with a factor A~ 4 and defining the dimensionless density po — X~ A pQ. The coordinate transformations and the 
identification relating J5J) with © are, 

r = p dx, 



= t+ I dx^^, (4) 



2 ' 
" 



2M po , 2 2 ^ 
J = — + -(c -v ). 

Notice that the coordinate r (or equivalently the dilaton <fi) from the point of view of the fluid plays the role of an 
extensive quantity. If the flux tube has a constant section it is proportional to the mass of the fluid contained between 
some reference initial point and point x. We stress the fact that the correspondence between 2D black holes and 
the acoustic metric © is exact. This is to be compared with what happens for 4D static black holes, for which the 
correspondence with the acoustic metric can be established only up to a conformal transformation of the metric Q . 

Eqs. allow to associate a one-dimensional slab geometry (and a corresponding acoustic black hole) to every 
2D black hole solution of the dilaton gravity model QJ. They can be also used to construct acoustic analogs of 4D 
spherically symmetric black hole solutions (e.g Schwarzschild). This is not surprising, being the causal structure of a 
4D spherically symmetric black hole completely encoded in its (r, t) sections. 
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III. 2D GRAVITATIONAL DYNAMICS VERSUS FLUID DYNAMICS 



In the previous section we have seen that given a generic 2D black hole we can always construct a 2D acoustic 
analog. However, the correspondence is purely kinematical. In order to formulate the correspondence at a dynamical 
level we need to compare the gravitational dynamics with the dynamics of the fluid, which determines po,v , c. 

At the dynamical level, the most striking difference between the two physical systems is that the gravitational 
dynamics has to be covariant under general coordinate transformations. In particular, this means that we have a 
huge redundancy of gauge degrees of freedom in the description of our 2D gravity system. These gauge degrees of 
freedom have to be fixed if we want to construct a dynamical acoustic analogue of the black hole. On the other hand, 
2D dilaton gravity is purely topological. We do not have any physical propagating gravitational degree of freedom. 
This leaves open the possibility that once the gauge freedom have been fixed, the dynamics for the resulting global 
gravitational degrees of freedom could match fluid dynamics. 

The field equations for the 2D dilaton gravity model Q are 

R = -X 2 ^, V M V„^ - \g^\ 2 V = 0. (5) 

We consider only static solutions and fix the diffeomorphisms invariance, choosing the Schwarzschild gauge for the 
2D metric 

ds 2 = -X{r)dr 2 + X- X {r)dr 2 . (6) 

Using Eq. © into the field equations J3J), one finds after some manipulations that the field equation for the dilaton 
4> are given simply by d<j)/dr = A, which implies that the dilaton is <fi = Xr. The field equations for the global 
gravitational degree of freedom X(r) become instead 

f.XV. (7) 
dr 

One can also easily check that the black hole J2J) is solution of Eq. Q . 

The fundamental equations describing fluid dynamics are the Euler and the continuity equations. The steady flow 
of the fluid giving rise to the one-dimensional slab geometry described in the previous section, is therefore governed 
by the equations: 

dvo dP dip _ . . . . . . . . 

Po«o-j 1- ~r + Po-r- = 0- po{x)Vo{x)A{x) = const., (8) 

dx dx dx 

where P is the pressure, ip denotes the potential for external driving forces acting on the fluid (including Newtonian 
gravity) and A(x) is the area of the section of the flux tube. We will consider separately the two cases: a) The external 
potential %p is non-homogeneous and the section of the flux tube is constant, b) The external potential is zero and 
the section of the flux tube does depend on x. 

Let us first focus on case a). Using the barotropicity condition P = P(po), passing from the coordinate x to the 
coordinate r as defined in Eq. 10} and defining the new variables 

X = — (C -V ), 

Y = p c, (9) 
F = In I — 



,Po, 

the Euler and continuity equations (JSJ take, respectively, the form 

dX _ ^dY x dF _ F dxl) 
dr dr dr dr 

±[Y(Y-X)}=0. (10) 

One can now easily realize that the fluid dynamics equations (|10|l can be made equivalent to the gravitational equations 
@ just by introducing the constraint 



(11) 



In fact, when constraint is enforced and the variable X is identified in terms of the gravitational variables as in 
Eq. © the Euler equation (fTTTf) matches exactly Eq. J7J). 

Simple counting of the degrees of freedom and of the fundamental equations of the dynamics reveals that we can 
consistently map the fluid dynamics into the gravitational dynamics imposing constraints. After gauge fixing, the 
gravitational system has only one degree of freedom, the function X parametrizing the metric in Eq. JBJ and one 
equation. The scalar field <fi is constrained to be proportional to r and therefore represents just the spacelike coordinate 
of the 2D spacetime. On the other hand, if the external potential is given, the fluid dynamics is characterized by three 
unknown functions c,po,vo and two equations, so that we have the freedom of imposing an additional constraint. If 
the external potential ip has been fixed, imposing the constraint is implicitly equivalent to an equation of state 
P = P(po) for the fluid. From this point of view, the equation of state for the fluid may be seen as a consequence of 
the gravitational dynamics. This is very similar to Jacobson's interpretation of the Einstein equation as an equation 
of state ^3 . Alternatively, if one first chooses a given equation of state the constraint (|llf) becomes an equation that 
determines the external potential ip. 

Two-dimensional gravity allows us to define the mass M of the black hole solution J5J as a scalar, which on-shell 
becomes constant 201 



M=i ( A 2 / r ( ,.),/o-(V,.)- 



(12) 



Owing to the correspondence between the dynamics of 2D dilaton gravity and the previously discussed constrained 
fluid dynamics, we also expect the existence of a function T of the fluid parameters X, Y, F, ip which is constant 
(independent of the flux tube coordinate x) by virtue of Eqs. 1)10(1 and (| 1 1 ft - This function is given by 



T = a 



[ {-XdF + 2e~ F dip) + 2Y - X 

J rn 



(13) 



where a and ro are arbitrary constants. One can easily check that dT/dr — follows from Eq. (|10J) and 1)11(1 . 

Eq. 1(13(1 allows us to find the acoustic counterpart of the black hole mass M. Fixing the constant a = A/2, using 
Eq. (jTT|) and reading from Eq. (0J the solution for X, X = J — 2M /A, one readily finds M = T. 



IV. BLACK HOLE THERMODYNAMICS AND FLUID DYNAMICS 



It is a well known fact that black hole thermodynamics follows directly from the gravitational field equations, once 
one has defined appropriately the thermodynamical parameters. In the case of 2D gravity the thermodynamical state 
of the black hole is completely defined by temperature T, energy (mass) M and entropy S. The mass is given by 
Eq. ((12(1 . The temperature can be defined in terms of the periodicity of the euclidean metric at the horizon, whereas 
the entropy can be computed using Noether charge techniques |2l| and turns out to be proportional to the dilaton 
evaluated at the black hole horizon. We have 

T=-^V(4> h ), S = 2<7r4> h . (14) 

47T 

The first principle of thermodynamics follows directly from Eq. Q). In fact, using Eqs. J5J one has S = 2irJ~ 1 (2M/ A), 
which in turn implies dM = (A/ 'Aii)(dJ '(r^) / 'd<ph)dS . On the other hand Eq. J7J implies dJ/dcp = V, which inserted 
in the previous equation and using Eq. ((141) gives dM — TdS. 

The correspondence between gravitational and fluid dynamics described in the previous section allows us to define 
straightforwardly the thermodynamical parameters, T a , M a , S a associated with the acoustic black hole. Moreover, the 
first principle of thermodynamics will become consequence of the fluid dynamics equations ((10(l . ((ll| ). 

Although the phenomenon of Hawking radiation for acoustic black holes is a rather involved issue p2l |. one can 
still define also for the acoustic black hole a temperature in terms of the " surface gravity" or the periodicity of the 
Euclidean section at the acoustic horizon c — vq. We have 



47r \ dr /„_, 47r 



(15) 



The above defined temperature coincides with the usual definition of temperature for acoustic black holes Q . Owing 
to the constraint (|1 1|) . in our 2D case we can also write down an other, completely equivalent, expression for the 
temperature, 
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For generic acoustic black holes it is not so easy to find a natural definition of mass. In the case under consideration 
we have a natural candidate for M a , given by the quantity (|13[) . which remains constant along the flux tube. Because 
it is a constant we can as well define M a as Eq. (| !!-{[) evaluated at the acoustic horizon , 



dr 



dr 



2e-^) +2y„ 7i : 



(17) 



where ro is chosen such that M a = 
gravitational black hole, S a = 2'K<j)} l 



M. The entropy of the acoustic black hole can be defined by analogy with the 
= 2-KTh- Using the coordinate transformation Q we get, 



S a = 2n 



dxp (x), 



(18) 



where we have fixed the arbitrary integration constant in such way that the integration domain is the whole range of 
the coordinate x outside the acoustic horizon. The entropy for the acoustic black hole is proportional to the total mass 
of the fluid outside the acoustic horizon. Notice that this definition corresponds to the intuitive notion of entropy as 
an extensive quantity, which counts elementary degrees of freedom. 

One can easily show, differentiating Eq. (|17|) and using Eqs. {Tfj)| and l|18[) that the thermodynamical parameters 
satisfy the first principle dM a = T a dS a - 

We can also show that the first principle can be derived directly from our constrained fluid dynamics. Evaluating 
the constraint (|llf> at the acoustic horizon (r = r^) we get XdJ(rh) = (2dY/dr — XdF/dr + 2e~ F dtp / dr) rh drh- Using 
Eqs. lUU) and (JTHJ together with dJ(r h ) = (2/X)dM, we find dM a = T a dS a . 



V. FLUX TUBE WITH NON-CONSTANT SECTION 



Until now we have considered a flux tube of constant section with an non-homogeneous external potential. The 
discussion can be easily extended to a flux tube with non-constant section. In this latter case the presence of a 
external potential is an unnecessary complication. For this reason we only consider the dynamics of a fluid with zero 
external potential and a non-constant section. In this case the fluid dynamics equations ©, written in terms of the 
variables defined in Eq. © are 



Introducing the constraint 



dr 



dr 



XV, 



(19) 



(20) 



Eqs. Ijl9(l can be made equivalent to the gravitational equations JJJ). 

The function that is constant owing to Eqs. I|19|) . (|2L)f> and can be therefore interpreted as the black hole mass is 



7 = a 



( (~XdF) + 2Y -X 

J rn 



(21) 



Also in this case, once A(x) has been fixed, the constraint H2()|l is equivalent to an equation of state for the fluid. On 
the other hand if the equation of state of the fluid is given, Eq. I|2(J|) determines A(x). In this latter case, Eq. I|2(J|) 
represents a "geometrical constraint" for the fluid flow of the same kind of that discussed in Ref. in relation to 
the realizability of acoustic horizons. 

The mass and temperature of the acoustic black hole are now 



A4 = I 



'•(-) 



dr 



2Y(r h 



T a = 



4?r/3o V dx dx J 



(22) 



whereas the entropy is always given by Eq. (|18fl . One can easily check the validity of the first principle of thermody- 
namics and that this principle is a consequence of the equations of the the fluid dynamics lj*H?)l . (|2*U|l . 
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VI. SOLUTIONS OF THE CONSTRAINED FLUID DYNAMICS 



Apart from its intrinsic theoretical interest, our approach could also have a direct experimental realization. Given 
a four-dimensional, spherically symmetric, gravitational black hole, one can derive by dimensional reduction the 
corresponding 2D dilaton gravity model QJ and related 2D black hole @- The 2D black hole solution is completely 
characterized by the function J and the mass M. Once J(r) and M are given, we can construct explicitly the 
corresponding 2D acoustic black hole and, at least in principle, realize the experimental set-up, which will produce 
the acoustic black hole. 

In order to have a well-defined physical system, we need to solve the fluid dynamics equations, i.e we need to know 
PojC, vq as function of the flux tube coordinate x. Equations @ and the second equation in i|10|) (or in i|19[l ) can be 
easily integrated 



X = J-™, Y= l -(X + ^X* + a 2 ) , (23) 

where a is a constant in case a), whereas a = const / A(x) in case b). 

To solve the dynamics completely, the constraint l|llfl (or I20H has to be enforced. This can be done in two 
nonequivalent ways. 1. One can choose a given barotropic equation of state for the fluid and the constraint (IJ2JJ) 
determines the external potential tp (the profile A(x) for the flux tube). 2. The external potential ip (the flux tube 
profile A(x)) is given and the constraint (|f I \ ( 1201) 1 determines the barotropic equation of state for the fluid. Because 
one usually operates with a fluid with a given equation of state, the physically relevant situation is obviously the first 
one. Therefore we will consider here only possibility 1.. To handle with a definite and common situation, we will 
consider a fluid with a constant speed of sound c. The barotropic equation of state is P = c 2 po, which expressed in 
terms of the variables 10 gives 

Ye F = c 2 . (24) 

Let us first consider the case of flux tube of constant section (case a)). Inserting Eqs. 124|) and l|23|l into the 
constraint (|If (I and using Eqs. l(9"jl. (|23|) . I|24|) one obtains the solution for vo(X), po(X) and ?/>(A), 

«o = ^, Po = ~Y, ^)=c 2 (|_-lny), (25 ) 

where Y(X) is given as in Eq. i|23|) . The acoustic horizon is located at X = (Y = a/2). At the horizon the density 
remains finite, po — a /2c. This is the best we can do in the general case. The calculations can be completed once a 
specific 2D gravitational model ( a specific function J) is given. Taking into account that po = dr/dx, the previous 
equation for p becomes a differential equation, which determines the coordinate transformation r = r(x). 

Let us now consider the case of a flux tube on non-constant section (our case b)). We can use Eqs. ill': ill to obtain 

A 2 cx — - -. (26) 

Y(Y-X) v ' 

Notice that Y > X. In order to determine Y(X) we use Eq. <|24[1 in the constraint 1|20|) to get 

dY XdY 

2 + = I. (27) 

dX Y dX V ' 

The solution of this equation is given in implicit form by 

X = 2Y In y, (28) 

where for simplicity we have set to zero the integration constant. We can now write the solution of our fluid dynamics 
equation in the form 



X Y 

«o = cWl-— , Po = — , (29) 

where Y(X) is implicitly defined by Eq. I|28|l . The horizon is located at X = 0, the subsonic region for X > 0. From 
very general arguments (see for instance Ref. 0) one expects that the acoustic horizon forms at a minimum of the 
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area of the flux tube section ( the so called Laval nozzle). Let us show that this is also true in our case. From Eqs. 
IJHJ and 1(23 one finds 

(3») 

where the prime denotes derivation with respect to X. On the other hand, Eq. I|27() tells us that at the horizon we 
have Y' = 1/2 and that in the subsonic (supersonic) region Y' < 1/2 (Y 1 > 1/2). It follows immediately that A' = 
at the horizon and that A' > (A' < 0) in the subsonic (supersonic) region, i.e that in the horizon region the flux 
tube must take the form of a Laval nozzle. 



VII. TWO EXAMPLES: 2D ANTI-DE SITTER AND 4D SCHWARZSCHILD BLACK HOLE 

In this section we will consider two part icular cases as illustration of our general procedure . Our first example is 
the 2D anti-de Sitter (AdS) black hole [23. The corresponding 2D dilaton gravity model Q has a constant potential 
V(<j>) = 2, from which follows J(r) = X 2 r 2 . The black hole solution |J2J) describes a 2D AdS spacetime with an horizon 
at r — \j2Mj A 3 . The corresponding acoustic black hole can be easily constructed using our general procedure. Let 
us consider separately the two cases of constant and non-constant flux tube section. 



A. Flux tube with constant section 



Considering a flux tube in presence of a non-homogeneous external potential and with constant section, the acoustic 
analog of the 2D AdS black hole can be easily found inserting X = X 2 r 2 — 2M/X into Eqs. Q25JI. Using the expression 
for po given in equation l|25|) one can also find the relation between the radial coordinate r of the gravitational black 
hole and the spatial coordinate x of the acoustic black hole. The < r < interior of the AdS black hole is described 
by the supersonic region of the acoustic black hole. At r = the fluid parameters po, vo and the external potential are 
finite. For r > v decreases and po grows monotonically. Also the external potential if) grows for r > but reaches 
its maximum at the horizon r = . At the acoustic horizon, vq = c and both the density po and the external potential 
remain finite. The asymptotic region r — > 00 of the gravitational black hole (the timelike conformal boundary of the 
2D AdS spacetime) corresponds to the subsonic region of the acoustic black hole. Going in this region vq — > , the 
fluid density p Q and fluid pressure P diverge as expected owing to the continuity equation. Conversely, the potential 
tjj of the external force acting on the fluid diverges logarithmically as — lnr. The behavior of the parameters if), i»o, Po, 
normalized to their horizon values and for a = 1 are shown if Fig. JlJ as a function of the variable X — \ 2 r 2 — 2M/X. 




X 



FIG. 1: Behavior of the normalized parameters ^/c 2 (normal line), «o/wo(r(,)(bold line) and po/pa{rh) (dashed line) as a 
function of the variable X for a flux tube with constant section and non-homogeneous external potential characterized by 
a = 1. The velocity and fluid density are normalized to their horizon values. The acoustic horizon is located at X = 



B. Flux tube with non-constant section 



If we use a flux tube of non-constant section to mimic the 2D AdS black hole, we get the typical shape of a 
converging/diverging Laval nozzle. We can find the solution for the fluid equations using X = X 2 r 2 — 2M/X into 
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Eqs. (|26|l . I|28|) l|29|) . In this case the fluid can be used to describe only the horizon region of the AdS black hole, 
but not the full spacetime. In fact, if A 2 of Eq. I|2t)|) has to stay positive X and because of Eq. I|28|) it follows that 
— 2/e < X < y/e and correspondingly 1/e < Y < \fe. The supersonic region of the acoustic black hole describes 
the region y/2MjX z — 2/(eA 2 ) = r m < r < rh = ^/2M/X 3 of the AdS black hole interior. Analogously, the subsonic 
region can be used to describe only the region y^2M/ A 3 — < r < Vm = -\/2M/A 3 + \/e/X 2 of the AdS black hole. 

The profile of the flux tube is that of the Laval nozzle depicted in Fig. @ , which also shows the behavior of the 
fluid velocity and fluid density, normalized to their horizon values, as a function of the variable Y defined in Eq. I|28|l . 
The flux tube section A starts at a finite value at r — r m decreases in the supersonic region of the acoustic black 
hole till it reaches its minimum Afr^) at the horizon. In the subsonic region < r < tm-i A grows and diverges at 
t = I'm- The fluid velocity decreases monotonically starting from vo(r m ) — \/3c in the supersonic region to reach first 
its horizon value vo(rh) — c and then vo(tm) = 0, as expected owing to the equation of continuity. Finally, the fluid 
density po, grows monotonically from po(r m ) = (ec) _1 in the supersonic region to the horizon value po(rh) = c _1 to 
reach po( r A/) = \Zec~ 1 in the subsonic region. 




FIG. 2: Behavior of the normalized parameters A/A(rh) (normal line), vo/vo(rh)(bo\d line) and po/po{rh) (dashed line) as a 
function of the variable Y for a flux tube with non-constant section. The flux tube section, the velocity and fluid density are 
normalized to their horizon values. The acoustic horizon is located at Y = 1 



Our second example is the 4D Schwarzschild black hole. By means of a dimensional reduction, retaining only the 
radial modes, 4D Einstein gravity can be described by the 2D gravity model Q with dilaton potential V(4>) = l/\f2<p 
and with A~ 2 identified as the 4D Newton constant G. The 2D black hole solution (J2J becomes, 

2M\ 




ds 2 = i^/2Xr - —j dr 2 + ( V2Ar - ^ ) dr z . ! TI ) 

This solution describes the 2D sections of the 4D Schwarzschild solution ds 2 = —(1 — 2MG/ R)dr 2 + (1 — 
2MG / R)~ 1 dR 2 . In fact the metric (|3I|I can be put in the Schwarzschild form by means of the Weyl rescaling of 
the 2D metric g^ v — > (l/^/24>)g flL , and the coordinate transformation r = (A/2)R 2 . Mass, temperature and entropy of 
the 2D black hole are invariant under Weyl transformations of the metric |24j . In fact the temperature and entropy of 
the 2D black hole l|31f) . given by T = X 2 /8ttM, S = (4ttM 2 /A 2 ), match exactly those of the 4D Schwarzschild black 
hole after setting A~ 2 = G. 

The acoustic analogue of the Schwarzschild black hole can be constructed, working with a flux tube of constant 
section putting X = V2Xr-2M/\ = XR- 2M/X into Eqs. $2~fy. (or into Eqs. ffity. $Zfy . |22J) if working with a flux 
tube of non constant section). 

The features of the acoustic Schwarzschild black hole are qualitatively similar to those of the previously discussed 
AdS acoustic black hole. 



VIII. SUMMARY AND OUTLOOK 



In this paper we have developed a general method for constructing acoustic analogs of 2D black hole solutions. The 
procedure can be also used for 4D black holes that admit an effective 2D description. In particular this is the case 
of spherically symmetric 4D black holes. We have also shown that after fixing the gauge degrees of freedom, the 2D 
gravitational dynamics is equivalent to a constrained fluid dynamics. The correspondence between gravitational and 
fluid dynamics has allowed a natural definition of the thermodynamical parameters mass, temperature and entropy 
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of the acoustic black hole. Moreover, we have shown that the first principle of thermodynamics follows from the 
constrained fluid dynamics. 

The results presented in this paper represent just a first step for finding condensed matter systems that mimic 
gravitational dynamics. It is not hard to identify the two main (and related) limitations of our approach. First, our 
discussion holds for static gravitational configurations. In many interesting situations, for instance the computation 
of the back-reaction for an evaporating black hole, one has to deal with non static solutions. The generalization of 
our method to the non static case may be extremely non trivial, in particular the fixing of the gravitational gauge 
degrees of freedom necessary to reproduce the fluid dynamics. The simplicity of 2D approach can be very useful also 
for treating the non static case. Second, in this paper we have considered a 2D gravitational model without matter. 
The absence of sources for the gravitational field represent a strong simplification for the dynamical equations. On 
the other hand, we need the presence of sources if we want to describe a realistic situation. Finding a fluid dynamic 
counterpart of the matter stress-energy tensor appearing in the gravitational field equations may not be a easy task. 

From a more technical point of view, if we work with a fluid with a given equation of state the practical realization 
of a physical system that mimics the 2D black hole requires a particular careful choice of an external potential (flux 
tube of constant section) or of the flux tube profile (flux tube of nonconstant section). 
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